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ABSTRACT 



GRS 1915+105 harbors one of the most massive known stellar black holes 
in the Galaxy. In May 2007, we observed GRS 1915+105 for ~ 117 ksec in 
the low/hard state using Suzaku. We collected and analyzed the data with the 
HXD/PIN and XIS cameras spanning the energy range from 2.3—55 keV. Fits 
to the spectra with simple models reveal strong disk reflection through an Fe K 
emission line and a Compton back-scattering hump. We report constraints on 
the spin parameter of the black hole in GRS 1915+105 using relativistic disk 
reflection models. The model for the soft X-ray spectrum (i.e.< 10 keV) suggests 
a = 0.56to;o2 and excludes zero spin at the 4(T level of confidence. The model for 
the full broadband spectrum suggests that the spin may be higher, a = 0.98]'^q;qJ 
(la confidence), and again excludes zero spin at the 2a level of confidence. We 
discuss these results in the context of other spin constraints and inner disk studies 
in GRS 1915+105. 



1. Introduction 

GRS 1915+105 was discovered by the WATCH instrument on board the Granat satellite 
on 1992 August 15 (Castro-Tirado et al. 1992) . It has since been classified as a microquasar. 
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a Galactic jet source with properties similar to quasars, but on a stellar scale. The mass of 
the central compact object is estimated to be 14±4 Mq (Grcincr ct al. 2001). The authors 
reported results on the optical counterpart that implied a K-M III donating star of ~ 1.2 
Mq, making GRS 1915+105 a low-mass X-ray binary. GRS 1915+105 is well-known for its 
extreme variability across all bands of the electromagnetic spectrum (Belloni et al. 2000). 

In 1998, BeppoSAX observations of GRS 1915+105 showed the clear presence of a 
broad emission line at 6.4 keV (Martocchia et al. 2002), i.e. the Fe Ka line. The Fe Ka 
emission line is a doublet consisting of the Kal and Kq;2 lines at 6.404 keV and 6.391 keV. 
Such emission lines arc likely produced by irradiation of the accretion disk by an external 
source of hard X-rays (an external X-ray source is required because accretion disks do not 
efficiently self-irradiate). It is generally thought that hard X-rays originate in a hot, diffuse 
corona that inverse-Compton scatters soft X-ray photons emitted by the accretion disk. 
Irradiation of the dense disk material by the hard X-rays then gives rise to a characteristic 
"reflection" spectrum that results from Compton scattering and photoelectric absorption 
(George & Fabian 1991). Photoelectric absorption within the disk leads to the formation of 
the fluorescent iron emission line at about 6.4 keV. In stellar-mass black holes, highly ionized 
iron lines are typically observed (Miller 2007). Because of its large cosmic abundance and 
high fluorescent yield, the Fe Ka line is the most prominent line in the X-ray reflection 
spectrum. 

Disk reflection is an important diagnostic of black hole spacctimc geometry and general 
relativistic (GR) effects (e.g. Fabian 1989; Laor 1991; Brcnncman & Reynolds 2006). It is 
believed that the iron line originates in the innermost parts of the accretion disk and thus is 
highly distorted by Newtonian, special relativistic beaming, and GR effects (e.g. Fabian, et 
al. 2000). This skewing explains the broad, skewed iron hne profile that is observed in some 
stellar-mass X-ray binaries, such as Cygnus X-1 and GRS 1915+105 (see Miller 2007). 

In our modeling of the X-ray spectrum, we make the assumption that the accretion 
disk extends down to the innermost stable circular orbit (ISCO) in the low/hard state. This 
may seem contrary to some models, in which the low/hard state is characterized by an 
inner disk boundary that is outside the ISCO. However, several pieces of evidence suggest 
that our assumption is likely to be accurate for GRS 1915+105. As follows, we observed 
GRS 1915+105 in the the low/hard state, which is characterized by an energy spectrum 
that is dominated by a power-law-like component. While there is a typical photon index 
of ~ 1.6 and a high-energy cutoff at ~ 100 kcV (Fender & Belloni 2004), there is also a 
very weak soft X-ray component (most likely associated with a thermal disk) (Miller et al. 
2006a). According to Nowak (1995), the total X-ray luminosity for systems in the low/hard 
state where distance and mass estimates are known is generally below 10% of the Eddington 
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luminosity. Esin et al. (1997), gave an archetypal model of how accretion flows change with 
state. In that model, the inner disk in the low/hard state is radially truncated at low mass 
accretion rates (corresponding to Lx/Ledd ^ 0.008) and an advcction-dominated accretion 
flow is present. However, the luminosity observed in GRS 1915+105 (L ~ 0.30 Ledd) is 
higher than what is theoretically calculated for a truncated disk in the low/hard state (Esin 
et al. 2007). This higher luminosity, and therefore higher mass accretion rate, is indicative of 
inner disk extending closer to the black hole (e.g. Esin et al. 2007: Fig. 1) and supports the 
notion that radially recessed disks in the low hard state may not always be the case. Recent 
analyses of GX 339-4, Cygnus X-1, SWIFT J1753.5-0127, and XTE J1817-330 have also 
revealed cool accretion disks that appear to remain close to the black hole at low accretion 
rates (Miller et al. 2006a, b, Rykoff et al. 2007, Tomsick et al. 2008). 

The radius of the ISCO is determined by the spin of the black hole: risco — 6.0 rg for 
a—0, and risco = 1-25 rg for an astrophysical maximum of a = 0.998 (where rg = GM/c^ 
and a, = cJ/GM^, see Thorne 1974, and Bardeen, Press, & Teukolsky 1972). The location 
of the ISCO determines the "strength" of the gravitational effects on the Fc K emission line 
(and broadband reflection spectrum), thereby influencing how broad and skewed the line 
proflle becomes. Therefore, determining the ISCO through line fitting serves as an indirect 
measure of the black hole spin parameter. 

Spectral fits to thermal emission from the accretion disk in black hole binaries provides 
an independent way to constrain black hole spin. In order to turn the flux observed from 
the disk continuum into a radius (and therefore a spin parameter), the source distance, 
inclination, and mass must be known. In addition, factors like the line of sight absorption, 
hardening due to scattering in the disk atmosphere, and the form of the hard spectral 
component can be important. In order to obtain reliable spin estimates using the continuum 
method, it is crucial for the source to have a strong thermal component (e.g. McChntock 
& Remillard 2009). For the case of GRS 1915+105, this thermal component is most easily 
observed while in the high/soft state. However, prior spin results that have been reported 
for GRS 1915+105 in this state are not consistent. McClintock et al. (2006) and Middleton 
et al. (2006) both used models for thermal continuum emission from the accretion disk 
to calculate values of d >0.98 and 0.721q;q2, respectively. In this paper, we present an 
independent analysis using the relativistic disk reflection spectrum to calculate the black 
hole spin in GRS 1915+105 in the low/hard state with the advanced cameras on Suzaku. 
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2. Data Reduction 

Suzaku observed GRS 1915+105 on 2007 May 7 starting at 14:40:22 (TT). The obser- 
vation duration was approximately 117 ks. The XIS pointing position was used. In order to 
prevent photon pileup, the XIS cameras were operated using the 1/4 window mode using a 
1.0 s burst option. The XISl and HXD/PIN cameras were used in this analysis. Two other 
XIS units were turned off to preserve telemetry and the fourth unit was run in a special tim- 
ing mode (which has yet to be calibrated). XIS on-source times of approximately 51 ks were 
achieved. This resulted in a dead-time corrected net exposure of 25 ks for the XISl camera 
(in 3x3 editing mode). A net exposure time of 57 ks was achieved using the HXD/PIN. 

For the soft x-ray data, the 2.3—10 keV energy range was used in order to avoid cali- 
bration problems near the Si K edge and because there are few photons at low energy owing 
to high line of sight absorption (nn ~ 4.0 x lO^^cm^^). Using the clean event file (3x3 
mode data) with the latest calibration databases at the time from the XIS data (CALDB 
20070731: version 2.0.6.13), we extracted the source light curve and spectrum with xselect. 
An annulus centered on the source with an inner radius of 78" and an outer radius of 208" 
(75 and 200 in pixel units, respectively) was used for the source extraction region because 
the center of the image suffered pile- up. An annulus centered on the source was also used for 
the background extraction region with an inner radius of 208" and an outer radius of 271" 
(200 and 260 in pixel units, respectively). We manually corrected for extracted region areas 
that did not land on the XIS chip. The XIS redistribution matrix files (RMFs) and ancillary 
response files (ARFs) were created using the tools xisrmfgen and xissimarfgen available 
in the HEASOFT version 6.6.2 data reduction package. The 3x3 mode event file was used 
to specify good time intervals and the data was grouped to a minimum of 25 counts per bin 
using the FTOOL grppha. 

For the HXD data, the latest calibration databases (CALDB 20070710) were used and 
the reduction process began with a clean event file from the PIN detector (energy range 
10.0—70.0 keV), which was at XIS aimpoint. In the PIN data reduction process, the PIN 
spectrum was extracted and deadtime was corrected by using the pseudo-events files. After 
the non- X-ray (NXB) background spectrum was extracted, the exposure time of the back- 
ground spectrum was increased by a factor of 10 since the NXB event file was calculated with 
a count rate 10 times higher than the real background count rate to reduce statistical errors. 
The cosmic X-ray background was simulated, modeled and added to the NXB spectrum in 
XSPEC version 12.5.0 as instructed in the Suzaku ABC guide. 
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3. Data Analysis and Results 

Using the X-ray spectral fitting software package (XSPEC v. 12.5.0, Arnaud 1996), 
we initially fit the Suzaku XIS spectrum of GRS 1915+105 with a simple absorbed power 
law model: phabs* (powerlaw) , with the phabs model component accounting for Galactic 
photoelectric absorption. All parameters were allowed to vary. The 4.0—7.0 keV range was 
ignored when fitting the data to obtain an unbiased view of the Fe K range. The fit resulted 
in an equivalent hydrogen column, nn, of 4.511q;o6 (in units of lO^^cm"^) and a photon 
index, T, of 2.09l|];[J} 

{"X^/^ = 1402/1282). That ignored energy range was then restored 
when forming the data/model ratio as shown in Figure 1. An asymmetric, skewed line is 
revealed in the exact way predicted for relativistic disk lines. In a similar manner. Figure 
2 shows the XIS and HXD spectra fit with a simple model consisting of a broken power- 
law modified by interstellar absorption: phabs* (bknpow). A constant was allowed to fioat 
between the two data sets. The 4.0—7.0 keV and 15.0—45.0 keV ranges were ignored when 
fitting in order to properly model the continuum. Note that for the HXD spectra, the model 
was fit over the 12.0—15.0 keV and 45.0—55.0 keV energy ranges. The fit resulted in n-a — 
4.35lg:2?, Ti = 2.031^2, and Fa = 2.59^^^ with an energy break at 7.7l^| keV {x^/v = 
1565/1312). Again, the ignored energy ranges were restored when forming the data/model 
ratio. Note that the residuals near 12 keV can be attributed to calibration problems near 
the edge of the detector. 

XSPEC provides a number of models for spectral line fitting. Models of particular 
interest to us not only allow for analysis of the observed broad iron line in our data, but 
account for the special and GR effects that physically go into creating such an emission 
profile around a black hole. Since determining the black hole spin of GRS 1915+105 through 
analysis of the iron line is our main priority, models that allow spin to be a free parameter 
in our fits are useful. Brenneman & Reynolds (2006) have calculated models (kerrdisk and 
kerrconv) that describe line emission from an accretion disk and include black hole spin 
as a free parameter, thereby allowing us to formally constrain the angular momentum of 
the black hole and other physical parameters of the system. The kerrdisk model describes 
line emission from an accretion disk, while the kerrconv model allows one to convolve a 
refiection spectrum with the smeared line function. 

Although general relativity permits the spin parameter a to have any arbitrary value, 
black holes can in principle have spin parameters -1 < a < 1. However, for simplicity the 
kerrdisk model only considers a black hole that has a progradc spin relative to the accretion 
disk that spins up to the Thorne (1974) spin-equilibrium limit, i.e., < a < 0.998. 

In addition to the spin parameter, we can specify and/or constrain nine other physical 
parameters in the kerrdisk model: (1) rest frame energy of the line, (2) emissivity index 
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within a specified radius, (3) emissivity index at radii larger than a specified value, (4) break 
radius separating the inner and outer portions of the disk (in units of gravitational radii), 
(5) the disk inclination angle to the line of sight, (6) the inner radius of the disk in units 
of the radius of marginal stability (rms), (7) the outer radius under consideration, (8) the 
cosmological redshift of the source {z=0 for our source), and (9) the normalization (flux) 
of the line in photons cm^^ s^^. kerrconv has the same basic set of parameters with the 
exception that it does not require an input line energy or a flux normalization parameter 
since it uses a kerrdisk kernel to smear the entire spectrum with relativistic effects. Of 
these parameters, we let 1, 2, 3, 5, and 9 vary freely. We bounded the rest frame energy for 
the iron hne between 6.4 and 6.97 keV and the disk inclination between 55 ° and 75 ° based 
on the figure presented in Fender et al. (1999: Fig. 6) from radio jet observations. The 
emissivity indices were fixed to be equal to each other meaning that we only consider one 
index, making the break radius meaningless. Parameters 4, 6 and 7 were frozen at 6.0 rg, 1.0 
rms, and 400 rms, respectively using the "standard" values of the kerrdisk and kerrconv 
models. By doing this, we assume that the ionization of the disk material within the ISCO 
is too high to produce significant line emission as has been shown by 3D MHD simulations 
(Reynolds & Fabian 2008). 

As shown in Fig. 3, our model is defined as: phabs* (kerrdisk + kerrconv*pexriv). 
A constant was allowed to float between the XIS and HXD data sets. The model compo- 
nent pexriv is an exponentially cut-off power-law spectrum reflected from ionized material 
(Magdziarz & Zdziarski 1995). The free parameters are the power-law photon index F, the 
fold energy or cutoff energy (in keV), and the scahng factor for reflection. Our results are 
given in Table 1 with the model being fit in the 2.3—10.0 keV and 12.0—55.0 keV energy 
ranges. Note that the residual feature at the upper limit of the XIS detector is the result of 
noise. Errors were calculated using the steppar command in XSPEC, which affords control 
over how the space is searched. For all parameters, errors refer to the 68% confidence 
level (Ax^= 1) unless otherwise stated. The reduced of our best-fit is 2345/2224. We 
measure a spm value oia^0.98to fi. Figure 5 plots the dependence of on the black hole 
spin parameter. A spin of zero is excluded at the 2a level of confidence. 

The disk ionization parameter = Lx/nr^, where n is the hydrogen number density) 
has a large effect on the resulting spectrum. Ross et al. (1999) demonstrated the effect on 
refiection spectra models for a range of ionization values ( 30< ^ < 10^). The model with the 
highest ionization parameter, ^ = 10^ was the best refiector but had negligible iron spectral 
features due to the disk surface layer being fully ionized at great depth, with the iron fine 
(specifically, Fe XXVI) not becoming dominant until r ~ 8. However, when ^ is reduced to 
3 X 10'^, less than half of the iron is fully ionized at the disk's surface and Fe XXV becomes 
dominant at r ~ 1. Most importantly, at this value an iron emission line that is Compton- 
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broadened becomes visible (Ross et al. 1999: Fig 2). Similar to this ionization value is our 
value of 5000, which is the upper limit allowed by the pexriv model. This "high" value 
was chosen so as to obtain the maximum amount of broadening due to Comptonization and 
further understand the dynamics of the system. In our model, maximizing or minimizing the 
disk temperature (10® and 10^ Kelvin, respectively) changed by less than 5, resulting in 
similar reduced to the best-fit value {x^ /f ~ 1.05). Effects to the parameter values (most 
notably the spin parameter) were negligible with the exception of the reflection scaling factor, 
which changed by < 0.06. Note that a thermal disk component is not required for our model 
and, since the disk itself is hkely to be relatively cool and the high Galactic photoelectric 
absorption prevents its detection, we have set a low disk temperature (see, e.g.. Miller et al. 
2006, Rykoff et al. 2007). 

Given the complexity of the source and it's extreme variability, we had set the iron 
abundance equal to the solar value to provide a more simplistic approach to constraining 
the spin, our main parameter of interest. Using the upper limit on the ionization value 
wc tested the effects of altering the iron abundance in the pexriv model. Ross & Fabian 
(2005: Fig. 5) show that increasing the iron abundance relative to the solar value serves 
to lower the continuum while enhancing the Fe Ka emission lines. Additionally, Lee et al. 
(2002) measured an iron abundance, Fe/solar, > 1.0 for GRS 1915+105, conveying that the 
overabundance may be due to dynamics of supernova in the microquasar's history. However, 
Ueda et al. (2009) were more definitive in fitting a spectrum with better statistics using 
solar values for the abundances (including iron). As such, after testing the effects of varying 
abundance, our pexriv model did not require an enhanced abundance (e.g. tripling the iron 
abundance increased by ~ 100). 

Our broadband model uses pexriv as its main refiection component. One difficulty 
with this model is that it does not include the effects of Comptonization on photoelectric 
absorption edges in the outer, most highly ionized layers of the disk. Being statistically 
inferior in this regime, it may not be ideal to use this model for high values of ^, which 
is the reason our value of 5000 erg cm/s is the upper limit. In order to check how our 
spin constraints might depend on the reflection model, we replaced pexriv with a newer, 
more self-consistent reflection model: reflionx (see Fig. 4). The reflionx model is a 
revised version of reflion (Ross & Fabian 2005) that describes reflection from accretion disks 
in black hole systems where the blackbody emission is at too low an energy to affect the 
Fe Ka emission. It has been used previously to constrain the spin in other sources such 
as GX 339-4 (Reis et al. 2008). The parameters of the model are the iron abundance 
(set to solar), photon index of the illuminating power law, ionization parameter, ^, and 
the normalization. In fltting reflionx, the HXD spectrum was neglected because a cut-off 
power-law is required whereas reflionx assumes that a simple power-law spectrum irradiates 
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the disk. Moreover, ref lionx is a pure reflection model that requires a simple power-law to 
represent the continuum. Note that the Fe K emission line is self-consistently modeled within 
the ref lionx model. Also, in order to obtain better constraints within the ref lionx model, 
the inclination was bounded between at 55 ° and 75 °. The model was fit in the 2.3—10.0 keV 
energy range. Letting a be a free parameter, the phabs* (powerlaw+kerrconv*ref lionx) 
model resulted in a spin parameter of 0.56 with a value of approximately 2124 and a 
reduced ~ 1-01 for 2100 degrees of freedom (see Table 2). Based on these results from 
the ref lionx model, we can state that a — is preferred and that a spin of zero is 

excluded at the 4(T level of confidence (see Fig. 6). However, it is important to note that the 
parameter values obtained with the ref lionx model are poorly constrained, which requires 
the pexriv results to be viewed with extra caution (see Sect. 4). 

Finally, we note the possible presence of an absorption line at approximately 7.4 keV 
(see Figure 1). It is not clear that the feature is real. When the broadband spectrum is fit 
phenomenologically, a Gaussian model for the line suggests a (single-trial) significance of Aa. 
The feature is not significant after fitting the relativistic emission line and disk reflection 
continuum. If the line is real, it could plausibly be associated with Fe XXV or Fe XXVI, 
implying a wind with a blue-shift as high as 0.1c. When a strong line is frozen as part of 
the overall spectral model, even modest constraints on the spin parameter of the black hole 
cannot be obtained. However, a wind in a hard state would be inconsistent with an apparent 
anti-correlation between winds and jets in GRS 1915-1-105 (Miller et al. 2008a, Nielsen & 
Lee 2009). 



4. Discussion 

In this paper, we have presented the flrst results from Suzaku observations of GRS 
1915-1-105 in a low/hard state. We have observed the Fe K emission line in GRS 1915-1-105 
in unprecedented detail due to the spectral resolution and fast readout modes of the XIS 
cameras. Our broadband spectral model suggests a spin of a = 0.98lo;oi) though a value of 
zero is only excluded at the 2a level of confldence. A different model, fltted only to the soft 
X-ray spectrum, results in a spin of a = 0.561o q2 with a value of zero excluded at Aa level 
of confldence. This effort to measure the spin of the black hole in GRS 1915+105 follows 
other recent efforts to measure spin in stellar-mass black holes using relativistic iron lines 
(Miller et al. 2008b, Reis et al. 2008/2009, Miller et al. 2009) and the thermal disk continuum 
spectrum (e.g. Shafee et al. 2006). 

Owing to its extreme behavior and the possibility of a relation to the spin of the black 
hole, it is particularly important to understand the black hole in GRS 1915-1-105. Middleton 
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et al. (2006) and McClintock et al. (2006) used similar thermal emission continuum models 
based on different data selections from RXTE to estimate the spin of the black hole in GRS 
1915+105. Middleton et al. (2006) determined the spin of GRS 1915+105 from a set of 16 
RXTE observations (approximately spanning 1994 to 1996) that were 16 s each. The authors 
chose spectra where there was a slow change between spectral states, specifically variability 
classes /3 and A (Belloni et al. 2000), which is why short time-scale spectral binning was 
used. They required that Comptonization contribute less than 15 % of the bolometric flux 
(or rather that the disk contributed greater than 85 %). Also, the rms variability had to be < 
5 %. For the final analysis, there were a total of 34 disk dominated spectra for 6 observations. 
The XSPEC models used were for a multicolor disk blackbody (diskbb), relativistic effects 
(bhspec), and thermal Comptonization (thcomp) that result in an intermediate spin value 
of a = 0.721^2 

For comparison, McClintock et al. (2006) used a much larger sample of RXTE data 
than Middleton et al. (2006) with a total observation time of 89 ks for RXTE and 13 ks 
for ASCA. The 22 observations selected had to meet 3 criteria: (1) weak QPOs, (2) rms 
continuum power < 0.075 rms, and a (3) disk flux > 75 % of the total 2-20 keV emission. 
The authors used a relativistic disk model (kerrbb2: a combination of kerrbb and bhspec 
that included a spectral hardening factor), which would alternately be combined with a 
power-law, Comptonization (comptt), and a cut-off power-law model (cutoffpl). These 
last 3 models represented the nonthermal tail component of the emission. They argued that 
the spin (a > 0.98) and mass accretion rate parameters were relatively unaffected by the 
model for the nonthermal tail component. 

Both sets of authors have identified several reasons that might explain the different spin 
values that they derived. McClintock et al. (2006) note the existence of uncertainties in disk 
structure at high luminosities. These uncertainties could possibly affect the Middleton et al. 
(2006) results, which used solely high luminosity observations (L > 0.3 Ledd)- McClintock 
et al. (2006) also suggest that Middleton et al. (2006) missed crucial low luminosity obser- 
vations because of small data samples - the derived intermediate spin value may be due to 
averaging over a wide range of luminosities (including a super-Eddington luminosity of L ~ 
1.45 Ledd)- However, Middleton et al. (2006) suggest that the McChntock et al. (2006) 
low luminosity data (L < 0.3 Ledd) does not allow for a low temperature Comptonized 
component (kTe < 3 keV). 

In summary, McClintock et al. (2006) finds a rapidly spinning black hole in GRS 
1915+105 based on disk continuum fits from RXTE and ASCA data. Middleton et al. 
2006) used similar assumptions (i.e. d ~ 12 kpc and inchnation angle ~ 66° (Fender et 
al. 1999)) and models, but had a smaller sample of RXTE data of GRS 1915+105 and 
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got an intermediate value of a = 0.72lo;o2- The behavior of GRS 1915+105 is particularly 

complex and the states identified in this source are more numerous and nuanced than those 
identified in other black hole binaries. The differences in Middleton ct al. (2006) can be 
partly attributed to the state classifications (i.e. A, B, C) for GRS 1915+105 as opposed 
to the black hole binary categories used by McClintock ct al. (2006). The spin estimate 
reported by McClintock et al. (2006) was obtained in a state that resembles the "thermal 
state" (i.e. high/soft state or "thermal dominant" state). Here, our observation was obtained 
in the low/hard state (which is mostly hkely "state C"). 

Previous work has suggested that for non-spinning black holes, the amount of Fc K emis- 
sion from within the ISCO is negligible: inside the ISCO, the ionization fraction sufficiently 
increases close to the black hole making it unlikely that there is significant Fe K emission 
within the ISCO since the ionization fraction substantially increases (Reynolds & Begelman 
1997). Reynolds & Fabian (2008) find that systematic errors on the spin due to emission 
from within the ISCO become smaller for more rapidly rotating black holes. Additionally, 
Shafee et al. (2008) determined that magnetic coupling is unimportant across the ISCO for 
geometrically thin disks. The best current theoretical works suggests that systematic errors 
are unlikely to bias our results for GRS 1915+105 (Reynolds & Fabian 2008). 

An r~^ disk cmissivity relation is typically assumed for line emission from a standard 
thin accretion disk with an isotropic source of disk irradiation (see, e.g, Reynolds & Nowak 
2003). Freezing the cmissivity index at 3.0 in the pexriv and reflionx models resulted in 
worse fits with A^;^ > 200 for 1 degree of freedom (F- value >180, > 8.0a level of confidence). 
Our current best-fit value for the cmissivity index in GRS 1915+105 is of interest because 
it implies that the reflected emission has a weaker dependence on the radius (r^^ °) than is 
typical. When the cmissivity index is high, as in the case for many spinning active galactic 
nuclei, there is significant iron line flux below 4 keV from material closest to the black hole 
and therefore high gravitational redshift. This is not true when the cmissivity index is low. 
Thus, spin determination for GRS 1915+105 case is being driven more by the detailed shape 
of the body of the iron line rather than the overall line extent. This could explain much of 
the difference between the pexriv and reflionx results. Moreover, a stable warp at the 
inner disk (e.g. due to Lense-Thirring precession) is a possible explanation for our different 
cmissivity index value in that it would reflect more than expected at larger radii. A link 
between Fe K emission and the phase of low-frequency quasi-periodic oscillations (QPOs) 
in the low/hard state of GRS 1915+105 may provide evidence for such a warp (Miller & 
Homan 2005). 

Martocchia et al. (2006) made two XMM-Newton observations of GRS 1915+105 in 
2004 during a long "plateau" (lower flux) state. The data was analyzed with EPIC-pn in 
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timing and burst modes, respectively. For the timing mode, they detected a broad excess 
at the energy of the iron hne and fit the data from 2.0—10.0 keV with the disk reflection 
model diskline (Fabian et al. 1989). Martocchia et al. (2006) found an inner radius of 
the disk in excess of 300 Rg for an unabsorbcd flux of l.OiO.l (10~^ crg/s/cm^) in the 
2.0—10.0 keV energy range. Although such a disk would have the advantage of providing a 
geometric explanation for the major X-ray spectral transitions seen in Galactic black holes 
(Done & Gierlinski 2006), the high (unabsorbed) flux (~ 3.3 x 10~^ erg/s/cm^) and hence 
high luminosity (~ 0.30 Ledd) we calculated in the 0.5—100 keV range for GRS 1915+105 
make it unlikely that the disk would be recessed to such a high value (e.g. Esin et al. 1997, 
Frank, King & Raine 2002, also see Miller et al. 2006a). We note that the high flux of GRS 
1915+105 may have caused photon pile-up in the EPIC-pn timing mode, thereby obscuring 
the breadth of the line and affecting the Martocchia et al. (2006) results. 

Based on VLA observations of relativistic ejections from GRS 1915+105, 70 ° ± 2 is the 
typical value used for the inclination of this source (Mirabel & Rodriguez 1994). However, 
in our broadband model the inclination angle pegged at 55 °, which is close to the lowest 
plausible value for this black hole binary (we set the bounds between 55 ° and 75 °, Fender et 
al. 1999: Fig. 6). Allowing the inclination to vary over a greater range (20°— 80°) resulted 
in our broadband model finding the minimum at 40°, which we disregard on the basis 
that it is not consistent with previous radio jet observations (e.g. Mirabel & Rodriguez 1994, 
Fender et al. 1999). Given our assumption that the spin axis of the black hole is parallel to 
the jet axis, our best-fit spin and inclination results follow from self-consistent assumptions. 

As mentioned previously, the resulting parameter values in the ref lionx model are not 
well constrained. To better demonstrate the dependence between parameters, particularly for 
the spin parameter, figures 7 and 8 show the 68, 90 and 99 % contour plots for inchnation 
versus spin for both the reflionx and pexriv models, respectively. As shown for the 
ref lionx model, there is a large degree of variation between the inclination and spin values 
whereas for the pexriv model, the spin and inclination values are more restricted. Due to 
the poor constraints in the reflionx model, the pexriv model results must be regarded 
with more caution. 

Among stellar-mass black holes, GRS 1915+105 displays uniquely rich and complex 
phenomena, both in X-rays and at shorter wavelengths. Its X-ray spectrum is similarly 
complex, and whether using the disk continuum or the reflection spectrum to constrain the 
spin, data seletion and modeling nuances can be important. Currently, our broadband model 
(pexriv) gives a spin parameter that is in agreement with the value reported by McClintock 
et al. (2006) (a > 0.98), while the upper limit on our soft X-ray spectrum model (reflionx) 
spin result is lower than the Middleton et al. (2006) spin value. In this effort to constrain the 
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Energy (keV) 

Fig. 1. — Data/model ratio obtained when the XIS 1 Suzaku spectrum of GRS 1915+105 was fit with 
a simple power-law model: phabs* (powerlaw) , that included photoelectric absorption. The 4.0—7.0 keV 
region was ignored when fitting the model. An asymmetric skewed line profile is evident in the way predicted 
for relativistic lines. 

spin of GRS 1915+105, we have considered only a single observation in a single state, though 
disk reflection is sure to be important in many states. Future observations with Suzaku and 
XMM-Newton in different states may allow the spin to be determined more clearly using the 
disk reflection spectrum. The thermal disk continuum is not evident in our observation, but 
future efforts that make joint use of the direct and reffected disk spectrum to constrain the 
spin (e.g. Miller at al. 2009) may provide another way forward. 

Acknowledgements 

We gratefully acknowledge the anonymous referee for helpful comments that improved this 
work. We thank the Suzaku mission managers and staff for executing our TOO observation. 
We thank Koji Mukai for helpful conversations. J.M.M. acknowledges funding from NASA 
through the Suzaku guest investigator program. M.C.M. acknowledges funding from NSF 
grant AST0708424. E.M.C. gratefully acknowledges support provided by NASA through 
the Chandra Fellowship Program, grant number PF890052. R.C.R. thanks STFC for grant 
support. This work has made use of the tools and facihties available through HEASARC, 
operated for NASA by GSFC. 



Table 1. XSPEC model 



Model component 



Parameter 



Value 



phabs 
kerrdisk 



nH (1022) ^cm 
Fe line energy (keV) 
rb (rms) 

Emissivity Index 
a (cJ/GM2) 
Inclination (deg) 
Ri 
Ro 

Redshift 

Norm, (photons/cm^/s) 
Photon Index 
Fold E (keV) 
Rel. Reflection 
Fe/solar 

Cosine of Inc. Angle 
Tdisk (K) 

Disk Ionization (erg cm/s) 
Norm. (photons/keV/cm^/s) 



,-2\ 



pexriv 



4 15+0-06 

^•-'-'-'-0.06 

g_40+o.o6 
(6.0) 

98+°-°^ 

55.0+20 

(1.0) 

(400) 
(0) 

0.010^ 



.002 
.002 
1 qR+0.03 
l.»U_o,o3 

-1.0 
-1.0 
h0.02 
-0.02 



38: 
0.38: 
(1.0) 
(0.45) 
(300,000) 
(5000) 

o o-l-O.lO 
'^■°-0.10 



2345/2224 



Note. — These are the best— fit parameters found using the 
model phabs* (kerrdisk+kerrconv*pexriv) with the XIS and 
HXD data. All of the kerrdisk and kerrconv parameters were 
tied and constant was allowed to float between the XIS and HXD 
data sets. The errors quoted in Table 1 are la errors unless 
otherwise noted. The values in parentheses were fixed during the 
fitting. The calculated luminosity (0.5—100 keV) (10"^^erg/s) for 
a distance d - 12 kpc (Fender & Belloni 2004) is 5.7to\o- 
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Table 2. XSPEC reflionx model 



Model component 


Parameter 


Value 


phabs 


nH (10^2) (cm-2) 


5 64+° °^ 


powerlaw 


Gamma 


/I /I +0.01 

^■44_o.oi 




Normalization 


6.78l°:°i 


kerrconv 


Index 


^•-■--0.1 




rb (rms) 


(6.0) 




a (cJ/GM2) 


0.561°:°^ 




Inclination (deg) 


69.3l|° 




Ri 


(1.0) 




Ro 


(400) 


reflionx 


Fe/solar 


(1.0) 




Gamma 


2.44lHi 




Xi 


10011° 




Redshift 


(0) 




Normalization (10~^) 


^•'-'-0.3 


x7^ 




2124/2100 



Note. — These arc the best— fit parameters found us- 
ing the phabs* (powerlaw+kerrconv*reflionx) model 
with only XIS data. The errors quoted are la errors un- 
less otherwise noted. The values denoted in parentheses 
are fixed parameters. 
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5 10 20 50 

Energy (keV) 

Fig. 2. — Data/model ratio obtained when the XIS and HXD Suzaku spectra of GRS 1915+105 were fitted 
with a broken power law model: phabs*(bknpow), that included photoelectric absorption. For the HXD 

spectra, the model was fit over the 12.0—15.0 kcV and 45.0—55.0 kcV energy ranges. The 4.0—7.0 kcV and 
15.0—45.0 keV regions were ignored when fitting the model. The residuals near 12 keV can be attributed to 
calibration problems near the edge of the detector. The curvature at high energy is a clear signature of disk 
reflection. 



r 




Energy (keV) 

Fig. 3. — Best-fit spectrum found through modeling of the iron line at 6.4 keV. Model: phabs* (kerrdisk 
+ kerrconv*pexriv) . 
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Spin Parameter 



Fig. 5. — The plot above shows the change in the goodness-of-fit statistic as a function of the black hole spin 
parameter, a for the model const . *phabs* (kerrdisk + kerrconv*pexriv) . Using the XSPEC steppar 
command, 20 evenly-spaced values of a were frozen and all other parameters were allowed to float freely to 
find the best fit at that spin parameter. The dotted lines indicate confidence intervals. 
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Spin Parameter 

Fig. 6. — The plot above shows the change in the goodncss-of-fit statistic as a function of the black hole spin 
parameter, d for the model phabs* (powerlaw+kerrconv*ref lionx) . Using the XSPEC steppar command, 
30 evenly-spaced values of a were frozen and all other parameters were allowed to float freely to find the 
best fit at that spin parameter. The dotted lines indicate confidence intervals. Please note that only the XIS 
data was used in this model. 
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Spin Parameter 



Fig. 7. — Inclination versus spin contour plot for GRS 1915+105 using the pexriv model. The 68, 90 and 
99 per cent confidence range for two parameters of interest are shown in black, red and green, respectively. 
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Fig. 8. — Inclination versus spin contour plot for GRS 1915+105 using the ref lionx model. The 68, 90 and 
99 per cent confidence range for two parameters of interest are shown in black, red and green, respectively. 
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